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Abstract 

We introdnce recursive corona graphs as a model of small-world networks. 
We investigate analytically the critical characteristics of the model, including 
order and size, degree distribution, average path length, clustering coefficient, 
and the number of spanning trees, as well as Kirchhoff index. Furthermore, 
we study the spectra for the adjacency matrix and the Laplacian matrix for 
the model. We obtain explicit results for all the quantities of the recursive 
corona graphs, which are similar to those observed in real-life networks. 
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1. Introduction 

For decades we want to know what a graph looks like. We want to reveal 
the principles of the networks’ behaviour covered by their complex topology 
and dynamics. We want to learn about how the network structure evolves 
over time and how it affects the properties of dynamical processes on it. 
For the nature of decentrality of real networks, it is hard to observe the 
networks directly. Instead we observe them by taking snapshots of their 
network structure and content and keep updating them. Yet this gives little 
information about the future, since many of them keep growing over time. 
Thus it is desirable to set up models to £t real networks in both structure and 
functionality. We can use the models to mimic real-life networks. We also 
expect that the properties of the models can be proven rigorously, thus we can 
hnd the relations between topological and dynamical properties of networks. 
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Even if it is hard to give closed-formed expressions for some quantities, it 
would be nice to make them tractable for convenience of estimation. 

Among various network models, the ER graph proposed by Erdos and 
Renyi [l| is the earliest one. It generates a random graph by choosing a con¬ 
stant probability for joining every pair of vertices in the network. The model 
exhibit interesting statistical properties and has been well studied by many 
people. However the model lacks some important properties of real-world 
networks. For example many real-world networks exhibit small-world prop¬ 
erty (3,1^ 1^ with their diameters growing logarithmically in the number of 
vertices, while maintaining a high clustering coefficient. The Watts-Strogatz 
(WS) model is a typical graph model with the small-world effect. Never¬ 
theless, because of its randomness, many of its properties cannot be derived 
precisely, for example eigenvalues of the adjacency and Laplacian matrices. 
Thus deterministic models are often used to mimic complex networks j^, 
since their structural and spectral characteristics can be determined 
analytically. In addition to the small-world effect, another important fea¬ 
ture of a network is degree distribution. Many real-world networks exhibit 
a heavy-tail distribution while some networks have an exponential distribu- 
The famous preference attachment 


tion [8|, 19|, llO 


111 Il2l. Il3l. Il4| scheme 


successfully described the growing process of networks with heavy-tailed de¬ 
gree distribution. This work, like the WS model , leads to a network with 
an exponential degree distribution. 

Recently graph (matrix) products have been applied to modelling graphs 
with the same properties as real-life networked systems, such as the Cartesian 
product [lij, dot product [l^ and Kronecker product [13, [Isl, 19, 21|. A 

merit of such methods is that the graph/matrix products facilitate estimate 
of the properties of the generated graphs. 

In this paper we introduce a recursive way to generate small-world net¬ 
works with an exponential degree distribution, based on corona product of 
graphs. We obtain exact solutions to many structural properties of the net¬ 
works. Moreover, we derive all the eigenvalues for their adjacency matrix 
and Laplacian matrix, which are provided in a recursive way. Based on the 
obtained eigenvalues, we calculate the number of spanning trees, as well as 
the Kirchhoff index of the networks. 
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2. Graph Construction 


In this paper, we use corona product to generate a small-world network 
model. Literatures about the corona product and its related graphs are 

M. 


Let G = {V{G), E{G)) be the embedded graph 


partly established 

of a network. Suppose the graph is undirected and has vertex set V{G) = 
{1, 2,..., iV} and edge set E{G) C V(G) x V(G). We dehne the number of 
vertices N as the order of the graph, and the number of edges as the size of 
the graph, denoted as M = |i?(G)|. 

Given two graphs Gi and G 2 , their corona product Gi o G 2 is dehned as 
follows. 


Definition 2.1. Let Gi and G 2 be two graphs with disjoint vertex sets. Gi 
has Ni vertices and G 2 has N 2 vertices. Their corona product Gi o G 2 is a 
new graph which consists of one copy of Gi nd Ni copies 0 /G 2 . The i-th 
vertex of Gi is joint by a new edge with every vertex in the i-th copy of G 2 . 

In this paper we investigate the case where G 2 is the g-complete graph, 
thus we give the dehnition of the recursive corona graph. 

Definition 2.2. Let Kg be the q-complete graph (q ^ 2), then the gth gen¬ 
eration of recursive corona graph (RCG) Gq{g + 1) is defined as the corona 
product of the previous generation of RCG Gq{g) and Kq. More formally, 
Gq{g -|- 1) is defined as Gq{g -|- 1) = Gq{g) o Kq, g ^ 0, with the initial 
condition Gq(0) = Kq. 

Figure [T] illustrates the construction process for a particular network 

C,{g). 


3. Structural properties 

In this section we derive several important quantities of the RCG, show¬ 
ing that it is an appropriate model for the small-world complex networks. 
Thanks to the deterministic feature of Gq{g), we can give exact expressions 
for the properties of the graph. We will give the explicit result for its or¬ 
der, size, degree distribution, degree correlation, average distance, clustering 
coefficient, number of spanning tree and Kirchhoff index. 

We denote by N{g) and M{g) respectively the order and the size of Gq{g). 
Next we show how to derive these quantities. Assume that the number 
of vertices and the number of edges that get newly generated at step g is 
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Figure 1: Construction of the recursive corona graph, showing the hrst three generations 

0fC4(ff) 

denoted as Ly{g) and Le^q)- Then it is obvious that we have Lv{g) = 
q X N{g — 1) for g ^ 1, which leads to the result of N{g) along with the 
initial condition iV(0) = q. With respect to the size of the network, we have 
LE^g) = g{g ~ l)/2 ■ N{g — l) + g ■ N{g — l),g^ 1, and the initial condition 
M(0) = g(g-l)/2. 

Proposition 3.1. The order and size of the graph Cq{g) = {V{g),E{g)) are, 
respectively. 

N{g) = q{q + IY (1) 

and 

M{g)^tq({q+1)0*^-2) ( 2 ) 

The average degree is 5{g) = —2(g + 1)“^ + g + 1 which tends to g + 1 
for large g. Note that many real-life networks are sparse and their average 
degree tend to a constant value. 

3.1. Degree distribution 

The degree distribution V{5) for a network is a function indicating the 
fraction of vertices with degree 5 over all vertices. The degree distribution is 
a very important characteristic of a graph. It is essential to the analysis of 
many other structural properties. 
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The cumulative degree distribution [8j is defined as 


S'=5 

which is often used to analyse the degree distribution of a graph. The 
quantity gives the fraction of vertices whose degree 5' is greater than or 
equal to 5. In addition, networks whose degree distributions are exponential: 
V{5) e have also an exponential cumulative distribution with the same 
exponent: 


y=5 s'=s '' 

Next we investigate the degree distribution of Cq{g). We find that at time 
= 0, the network has q vertices of degree q — 1. Now we study the degree 
of some vertex v at step g. Let the value be 6y{g), we look in details about 
how the quantity evolves. We assume that vertex v is added to the network 
at step g^, {g^ ^ 0). For any g^ > 0, there are Sy{gy) = q, where q — 1 edges 
link to other vertices in Kg and the other 1 edge links to Cq{g — 1). At every 
step every existed vertex increase its degree by q. 

Theorem 3.2. The cumulative degree distribution of the graph Cq{g) follows 
an exponential distribution: Vcumid) (g + 1) 

Proof. The degree of vertex v at step g, denoted as Sy{g), can be written as 

Sy{g + 1) = 6y{g) + q . (4) 


Thus we have 

Svig) = q ■ {g - gy + l) ,igy > 0), and 

Sv{g) = q-{g + l)-l,{gv = 0). (5) 

This means that the numbers of vertices with the degree equal to q, 2q, ■ ■ ■ , gq^ {g+ 
l)q — 1 are, respectively, + 1)^“^, q^{q + 1)^“^, • • • ,q^,q- 
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For a certain value of degree 6, we have Pcum(<5) = (^-/V(O) + X]g'=i ^vig ')) /N{g) 
where 6 = \_g — {S + q)/q\. Therefore we can find 


Vcura{S) = 


q{q + iy 

1 


9 

g' = l 

9 

g 


-1 




For large g we have 


{q + iy + 

= (g + iy-^- 

Pcum(5) = (g + i)'""~(g + i)"^+'. 


( 6 ) 


□ 


3.2. Degree correlation 

One important parameter for the degree correlation is the average degree 
of adjacent vertices of all vertices which is referred to as any vertex with 
degree 6. We denote the parameter by fcnn(<5). If A;nn(^) increases with 6, 
this means that the vertices have a tendency to connect to vertices with a 
similar or larger degree. In this case we claim the graph to be assortative. 
The considered value of vertices with degree 6 which is generated at step g^ 
can be written as 


Qv-l 


\ 9p—^ 

9 \ 

+ Lv{gv){q - D)q{g - + l) + L^{g^) ^ q-q{g-i + D)\ 

i=gv+l / 

1 , , , l + q-2{l + q)-s^+^ 

—i;q{g — gv + D —jr • 

2 q{g - gv + D 


( 7 ) 


According to Eq.(l5]), we can express it as 


, ,,, 1, g + l-2(l + g)"^+^ 

kUS) = -{q + <5) + 2^-- 
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Since q ^ S ^ q{g + 1) — 1, therefore we have 0 < (1 + g)? ^ < 1, which 
means 

t.nW “ + ■?) + f (8) 

As for the initial vertices, everyone of them has the same distribution on the 
degrees of its neighbors. It leads to 


=^—~ + 1) - 1) + g - * + 1) 

{9 + 2)g {g + 2)g — 2 


(9) 


which yields 


2(gg + g- 1) ’ 


^nn(<^o) — 77 ( ^0 + O' + 


1 -g 


( 10 ) 


By checking the results we can see that the considered graph is assortative. 


3.3. Average distance 

Given a graph G = (V, E), its average distance or mean distance is dehned 
as: /i(G) = |y(G)|(iy(G)|_i) ■ Y.u,v(iV{G)d{u,v) where d{u,v) is the distance 
between the pair of vertices u and v. 


Theorem 3.3. The average distance of graph Cq{g) is 


KCq{g)) 


(g+1) g (2gg^(g +l)^g ^ + (g + 1)3 + (g - 2) ((g + 1)^)^) 

g(g +1)9-1 


( 11 ) 


Proof. To begin with, we assume that the summation of distances between 
all pairs of vertices in Cq{g) is D{Cq{g)). The sum of distances between all 
pairs (m, v) where u belongs to vertex set U and v belongs to a disjoint vertex 
set V, is denoted as D{U,V). 

In order to utilize the recursive construction process to the recursive 
corona graph we classify the vertex pairs in Cq{g) into 4 different categories 
W, X, Y and Z. The sum of the distances for the 4 categories is denoted as 
Sw, Sx, Sy and Sz, respectively. 

Category W refers to the pairs within the same Kq that we add to the 
network at step g. Category X refers to the pairs u,v where u is selected 
from one of the N{g — 1) KqS added at step g and v selected from any other 
Kq added to the network at the same step. Y refers to the pairs where u is a 
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new vertex and f is a vertex in Cq{g — 1). As for category Z, it indicates the 
pairs where both u and v are from the previous generation of graph Cq{g — 1). 
Thus we have the following equations: 


D{Cq{g)) 

Swig) 

Sxig) 


Svig) 


Szig) 


Swig) + Sxig)+ SYig) + Szig), (12) 

^h^.iV(g-l), (13) 

(dij + 2) ■ 

i<j 

q^Szig) + q^Nig - 1) • (iV(g - 1) - 1), (14) 

(djj + 1) ■ g 

2qSzig) + qNig - 1) • (iV(g - 1) - 1) + iV(g - 1) • g ,(15) 
DiCqig-1)). (16) 


Combining these recursive expression together we have: 


D(C,{g)) = (<, + l)’B(C,(g - 1)) + Ig'^ (2g(g + l)»-‘ - l) (g + 1)». (17) 


with the initial condition Zl(Cq(0)) = g(g — l)/2 we get the result 

0(C,(g)) = ig (2gg^(g + 1)^»-' + (g + 1)» + (g - 2) ((g + !)=)») . (18) 


Dividing D(Cg(g)) by i) yjgi^jg Eq. fITT]) . For large g, we have p(Cq(g)) ~ 

2 g 2logqNig), which increases logarithmically with the network order. 

□ 


3.4- Clustering coefficient 

Clustering coefficient is another crucial quantity used to characterize 
network structure. Many works about determining clustering coefficient and 
its related quantities are done on both graph models and graphs in reality j^, 


The clustering coefficient of vertex v is dehned as the following quantity 


civ) 


2ev 

6yiS^ - 1 ) ’ 

8 


(19) 








where is the number of edges between the neighbours of vertex v. The 
network clustering coefficient C{G) is dehned as the average of c{v) among 
all vertices. That is, 


c(G) 


1 


c(n). 

v&V(G) 


( 20 ) 


Theorem 3.4. Let v{6) be a vertex in Cq{g) whose degree is 6. Except the 
initial vertices, its clustering coefficient is 

c(n(5)) = |^^^. (21) 


Proof. Let us review the intermediate result in calculating the degree distri¬ 
bution that the number of vertices with degree q,2q, - ■ ■ , gq, {g + l)q — 1 are, 
respectively, -|- 1)^“^, q^{q -|- 1)^“^, • • • ,q‘^,q- Except the initial vertices, 
the clustering coefficient of other vertices follow the same rule, that is, a 
vertex with degree k ■ q has k ■ q neighbours, which are evenly distributed 
in k clusters. Each cluster forms a complete graph Kg. Thus the clustering 
coefficient of vertex v is derived as 


kq{q-l)/2 q-1 

= m-TTTTt = 1 -7) idv > 0) , 


kq{kq — l)/2 kq — 1 


and for initial vertices 


, ^ (g- l)(g- 2 ) +^g(g-1) , 

= (,(9 + 1 ) - !)(,(«+ 1 )- 2 )’ '«• = ■ 

Theorem 13.41 is naturally gained. 


( 22 ) 


(23) 

□ 


Theorem 3.5. The clustering coefficient of RCG network Cq{g) converge to 


c{Cq{g)) ~ 


g +1 


-<h 




(24) 


when the network order is high enough. In the expression, <h is the Lerch 
transcendent function. For large q, the clustering coefficient tends to 1. 
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Proof. From Eg. 020 p and 0231) we can obtain the clustering coefficient of 

CM 


ciCqia)) = 


1 


Nig) 
g-i 
q + l 


E 

.k=l 

V 


q-1 

kq — 1 




jq - l)igq + q - 2) 

(9(5 + 1) - 1)(9(9 + 1) - 2) 


1 ^ 


1 


q+l 


<?+l 




1 ^ gg + g - 1 ^ (g + 1) 


g +1 


g 


gg + g - 1 


(25) 


When g is large enough, we obtain 

q - 1 


ciCM) 


$ 


g + 1 \q + l 


1 ^ g-1 

q 


(26) 


which is high and tends to 1 when q is also large. Figure Ogives the clustering 
coefficient of some networks while they grow. □ 



Figure 2: Clustering coefficient of some recursive corona graphs. 
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3.5. Spanning trees 

Next we derive the number of spanning trees in graph Cq{g). 

Theorem 3.6. The number of spanning tree of Cq{g) is 

A'.r(C,(s)) = <;’-"(9+l)‘-'’(('i + l)'’'')'''^‘’' ■ (27) 

Proof. According to the Cayley’s formula j^, the number of spanning trees 
of a complete graph Kq is equal to Since all vertices of a Kq added to 
the graph are connected to a vertex in the original graph, these q+1 vertices 
consist of a new complete graph iCq+i- Therefore we have the following 
recursive relation of the spanning trees of Cq{g): 

NtriCqig)) = {{q + NUg - 1). (28) 

Together with the initial condition N^r = we can derive the expression 
for Ntr] 

Ni,{C,(g)) = r-H‘1 + 1)'“’ ((9 + !)»-■)''’+"" . (29) 

□ 


3.6. Kirchhoff index 

Resistance distance is an important character of a graph, which can imply 
many of its dynamic properties. The Kirchhoff index 28| of a graph refers 
to the sum of resistance between all vertex pairs in an associated electrical 
network obtained from the graph by replacing each edge of the graph by a 
unit resistance. Denote the effective resistance between vertices i and j as 
r{i,j) or Tjj, then the Kirchhoff index i?Kr of graph G is dehned as 


flK,(G)= 




(30) 


ijeViG) 

i<j 


We denote the Kirchhoff index of graph Cq{g) by RKr{Cq{g)). 

Theorem 3.7. The Kirchhoff index of Cq{g) is 

RKr{Cq{g)) = {q\2g + 1) - 2g - 1) (g + + q{q + 1)^-^. (31) 
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Proof. We denote by r{U, V) as the sum of all effective resistance between 
pairs (m, n) in which u and v belong to two disjoint vertex set U and V respec¬ 
tively. Similar to the method we used in calculating the average distance, we 
classify these pairs into 4 categories W, X, Y and Z, where the dehnition is 
exactly the same as used in calculating the average distance. Then the sum 
of the distances for the four categories is denoted as Rw, Rx, Ry and Rz- 
We have the following equations: 

RKriCqig)) = 

Rwig) = 

Rx {g) = 


Rrig) = 


Rz{g) = 

which yields 

RKr{Cq{g + 1)) = if^qi^q + 1)® “ 1) (^ + 1)^ + (O' + Rk-AC q{g)). (37) 

Notice that the effective resistance between vertices v and m in a complete 
graph Kq is 2/g since the potential between any other vertices are identical, 
if we impose potential difference between u and v. 

Along with the initial condition Rk{Cq{R)) = g — 1 we can deduce 

RKr{Cq{g)) = {q^-2q-l){q + lf3-^ + q{2q^g{q + iy + q + l){q + iy-^ 

= {q^-2q-l + 2gq^){q + 

= {qY2g + 1) - 2g - 1) (g + + g(g + l)^-'. (38) 

This completes the proof. □ 


Rw{g) + Rx{g) + Rvig) + Rz{g ), 

g(g - 1) 2 


g +1 


N(g-l), 


E (^M+2-^)-g 


i^3^y{Cq{g-l)) 
i<j 


q^Rz{g) + 


2 • q^Njg - 1) • {Njg - 1) - 1) 

g + 1 


2qRz{g) + q • ^-^N(g - 1) • {N{g - 1) - 1) 


+q 


g + 1 

RKr{Cq{g — 1 )) , 


N{g - 1), 


(32) 

(33) 


(34) 


(35) 

(36) 
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4. Spectral analysis 

By convention the (unweighted) adjacency matrix A{G) of a graph G is 
defined as a. N x N matrix with the entry ajj representing the number of 
edges incident with endpoints i, j. The degree matrix of G, D{G) is defined as 
a diagonal matrix with its ith entry on the main diagonal equal to the degree 
of vertex i. We call L{G) = D{G) — A{G) the Laplacian matrix of graph 
G. These matrices determine the structure of the graph, and the eigenvalues 
of A{G) and L{G) are sensitive to many of the structural properties, which 
have remarkable impact on the dynamic processes superimposed upon the 
network. 


Definition 4.1. Given A{Gq{g)) the adjacency matrix of Gq{g), we define 
the spectra of Gq{g) as 

^{Cqih)) ■= (^{g) = (Ai(^), X2{9) ■ ■ ■ Kig)) • (39) 

Similarly, we have: 

Definition 4.2. Given L{Gq{g)) as the Laplacian matrix of Gq{g), we define 
its Laplacian spectra 


S{Cq{g)) := S{g) = (7i(^), 72(^), • • • ,7n(^))- (40) 


4-1. Spectra of Adjacency Matrix 

Theorem 4.1. The relation between a{g) and a{g — 1) is 

1. i)+<? i±\/(g 1 >^iia 1)) ^ multiplicity 1 fori 

1)^“^ and 

2. —1 G a{g) with multiplicity {q — l)q{q + 1)^“^. 


The result is a corollary of results in [22|, [29 


1, ••• ,q{q+ 


4-2. Spectra of Laplacian Matrix 

Theorem 4.2. The relation between S{g) and S{g — 1) is 

1. g 1 fori = 

1, • • • , q{q + 1)^“^ and 

2. g + 1 G S{g) with extra multiplicity (g — l)g(g + 1)®“^. 
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Note that in the hrst part 7 ^( 5 ' — 1) = 0 will generate an eigenvalue equal 
to g + 1 with multiplicity 1 in iteration g. So the actual multiplicity of g + 1 
is (g — l)g(g + 1)^“^ + 1 for any g ^ 1. The proof of Theorem 14.21 is evident 
using methods in j^, 23, 2^- For convenience of the following discussion we 
give a similar proof here: 


Proof. The Laplacian matrix of Cq{g) is 


/ L(C,(g-l)) + g 4 


HCq{g)) = 


V 


-In 

-In 


-In \ 


{L{Kq) + /,) 


/ 


(41) 

Let hi, • • • , yAr(g_i) be the Laplacian eigenvectors of Cq{g — 1) correspond¬ 
ing to the eigenvalues 7 i(g - l), 72 (g - 1 ), • • • , 7 iv(g-i)( 5 ' - 1 ), respectively. 
For i = 1, • • • , iV(g — 1), let 

, _ li{g -l)-fg + l-fV {ii{g -l) + q + lY - - 1) 

2 ’ 

7 _li{g - 1 ) -f g + 1 - {ii{g -l) + q + lY - ^-ii{g - 1) 


Note that 0* are Laplacian eigenvalues of Cq{g) corresponding to the 
eigenvectors 


i Yi \ 

9{(l>i)Yi 

\ 9{(t>i)Yi ) 


( \ 

9{k)Yi 

V 9{.k)Yi ) 
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respectively. In fact 0* is obtained by solving 


/ Yi \ 

g{(t>i)Y^ 


L{C,{g)) 


( Yi \ 

9{4’i)Yi 


\ 9{4’i)Yi / 

/ L{Cg{g - 1)) + qin \ 
-In 

\ / 


\ 

Y^ 


9{(t>i)Yi ) 

( 

<? 

+ E 

i=i 


-In \ 

{L{Kq) + In)ljln 

{L{Kq) + In)qjln j 


9{(t^i)Yi 


/ 'Jiig- 1) +q-q-gi4>i) \ 

-l+g{(pi) 


\ -1 + 5(00 / 


Thus we can derive the following equations: 


(42) 


t(5 -l) + q-q- 5(00 = 0*, (43) 

-1 + 5(00 = 0i-5(00. (44) 

From Eq. (l44)) we can obtain that 5(00 ^ 0. Therefore we can substitute 
Eq. fHT)) into Eq. fH^ . we have: 


(705 - 1) + 5 - 00 (0i - 1) = -5, (45) 


which leads to the result of the hrst part of the theorem. 

If the Laplacian eigenvalues ui = 0, 1^2 = I's = ■■■ = I'q = q oi L{Kq) 
are correlated with the eigenvectors Zi,Z 2 ,-- - ,Zq, respectively, then for 
j = 2 , • • • , 5 , we have 

This completes the proof. □ 

Next we use the results of the Laplacian spectra to prove Theorem 13.61 
and 13.71 First we give an alternative proof of Theorem 13.61 
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Proof. It is known that the number of spanning tree of a graph G has the 


following form j30|, 


n N 

^=2 P 


Ntr{G) = 


N 


(47) 


where N is the number of vertices and r* refers to N eigenvalues of the graph 
G. Given the graph is connected, let ri be the unique zero eigenvalue, then 
Ti, i = 2, ■ ■ ■ ,N are N — 1 non-zero eigenvalues of the graph Laplacian. 

Theorem 14.21 tells that the Laplacian spectrum of Gq{g) consists of two 
parts. For the hrst part, we can derive from Eq. dTSj) that, in iteration g, 'ji 
in Gq{g — 1) generates two eigenvalues 0* and (pi which are subject to the 
relations 0*0* = 7 * and 0*4-0* = 7 *-|-g-|-l. In particular the trivial eigenvalue 
7 i = 0 generates 0* = g -|- 1 and 0* = 0. As for the second part, there is an 
eigenvalue 7 = g -|- 1 with multiplicity (g — l)g(g -|- l)^“hWe denote by S{g) 
the sum of all non-zero eigenvalues of L{Cq{g)) and by T(g) the product of 
all non-zero eigenvalues of L{Cq{g)). Then we can obtain 


Nig) Nig- 1 ) 

F(5) = n = (9 + ^ (pipi 

i =2 i =2 

Nig-i) (48) 

= (^ + l)(9-ik(9+iF- +1 7*(g-l) 

i =2 

= T{g - 1) • g(g + i)('?-i)9(<?+i)«-i+i , 

Eq. (H 8 ll and the initial condition T(0) = q‘^~^ yield 

T(g) = g'?-i(g + l)(9-b((9+iF-i)+9 . ( 49 ) 

Therefore 

Ntr{Gq{g)) = q^-\q + . 

□ 


The result is equivalent to what we derived using combinatorial method. 
In the following we give an alternative proof of Theorem 13.71 using the 
spectra information. 


Proof. The Kirchhoff index of a graph G can be expressed as [32, l33 


RKr{G)=N-^-, 
A-n P 


( 50 ) 
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where N and r, are the same as the previous definition. Let 


N N N N ^ 

^=En^'‘=n^*-E-- (51) 

j=2 i=2 i=2 '* 

i¥=j 

We can follow the clue of the previous analysis by separating the eigenvalues 
of its Laplacian matrix into two parts. Recall that the eigenvalues of the 
Laplacian consist of two parts and as defined by Theorem 14.21 As¬ 
sume that \ {0}. For the first part of the eigenvalues of L{Cq{g)), 

we denote them as (/)j and (j)i, i = 1,2, ■■■ , N{g — 1). Suppose the original 
eigenvalue in L{Cq{g — 1)), which is correlated with and (pi, is 7 *. Then 

A(^) = A(i)T(2)+ R(2)tW, (52) 

where 

n ’ (53) 

^(2) = JJ 7 i, (54) 

7i6r{2) 

E n (55) 

7jer(2) 

i¥=3 

and 

7^^^ = n ■ (55) 

7ier(i) 
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Accordingly we can obtain 

S^^^={q + 1) ^ 


^2 'ri 

i=2 \ j=2 




N{g-1) 

+ (ki^i 

i=2 


Nig-1) 


+ n 


/ 


i=2 


N{g-l) /N{g-l) \ 

= (^ + 1 ) n • (</'j + 

V s? 

Nig-l)N(g-l) 

=T{g-l) + {q+l) 'yi{g-l)-{'yj{g-l) + q+l) 

5? 

N{g-l)N(g-l) 

=T{g-l) + {q+l) Y n 7i(5'-1) •7i(5'-1) 

j=2 i=2 

N(g-l)N(g-l) 

+ iQ + '^) Y n 7i(5'-1) • (g +1) 

7=2 

= (1 + (<?+!)• {q{q + - 1)) • T{g - 1) + {q + 1)^ • S{g - 1). 

Y(2) = [q-\- ^ 

S(2) = (q- l)q(q + i)(q-l)9(g+l)«-i+s-2 ^ 


(57) 


(58) 

(59) 


= (5 + 1 )-T( 5 - 1 ). (60) 

Therefore we obtain a recursive relation of S{g). Considering the initial 
condition S'(O) = {q — we derive 

S{g) = q'^-Hq + l)(''-i)('?+i)«+9-.-i (((25 + 1)^3 _ 25 _ 1) (5 + ly + 5(5 + 1)) , 

(61) 

thus 

RKr{Cq{g)) = {q + iy~^ {{{2g + l)q^ - 2q - 1) {q + ly + q{q + 1)) 

= {q^{2g + l)-2q-l){q + iy3-2^q^q^yg-y ( 52 ) 

□ 

Note that Eq. fl6^ is consistent with Eq. (HT} . For large g, the Kirchhoff 
index displays the following leading behavior: 

RKr{Cq{g)) ~ g ■ N{gf ~ N{gf log,+i N{g) , (63) 
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5. Conclusion 

In this paper, we have introduced a deterministically growing model to 
generate small-world graph, by using the corona product. The advantage of 
such a model is that many of its properties can be solved exactly. We have 
derived explicitly many structural quantities of the small-world model. We 
have also found the eigenvalues for the adjacency matrix and the Laplacian 
matrix of the model. In future, the properties of various dynamical processes 
taking place on the small-world model deserve study. 
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